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Abstract 

Invariant manifolds facilitate the understanding of nonlinear stochas- 
tic dynamics. When an invariant manifold is represented approxi- 
mately by a graph for example, the whole stochastic dynamical system 
may be reduced or restricted to this manifold. This reduced system 
may provide valuable dynamical information for the original system. 
The authors have derived an invariant manifold reduction or restric- 
tion principle for systems of Stratonovich or Ito stochastic differential 
equations. 

Two concepts of invariance are considered for invariant manifolds. 

The first invariance concept is in the framework of cocycles — an 
invariant manifold being a random set. The dynamical reduction is 
achieved by investigating random center manifolds. 

The second invariance concept is in the sense of almost sure — an 
invariant manifold being a deterministic set which is not necessarily 
attracting. The restriction of the original stochastic system on this 
deterministic local invariant manifold is still a stochastic system but 
with reduced dimension. 

Key Words: Stochastic differential equations; invariant mani- 
folds; random center manifold reduction; almost sure invariance; method 
of characteristics. 
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1 Introduction 



Invariant manifolds provide geometric structures that describe dynamical 
behavior of nonlinear systems. Dynamical reductions to attracting invari- 
ant manifolds or dynamical restrictions to other (not necessarily attracting) 
invariant manifolds are often sought to gain understanding of nonlinear dy- 
namics. 

There have been recent works on invariant manifolds for stochastic or 
random ordinary differential equations by Carverhill [§] , Wanner j3U] , Arnold 
[3J, Boxler [HUH], and Mohammed [21], among others. These authors use the 
(sample-wise) cocycle property for the solution operator of the stochastic 
differential equations, the Osledets' multiplicative ergodic theorem [3J, and a 
less-physical but technically convenient random norm, to prove the existence 
of invariant manifolds. The construction of a random norm needs the knowl- 
edge of Oseledets spaces (a kind of eigenspace in random linear algebra) as 
well as Lyapunov exponents, both are hardly ever available; see [3J, P- 191 
and p. 379. Random norms are not realistic in this sense, and thus represen- 
tations of invariant manifolds and the dynamical reductions are difficult to 
achieve when random norms are used. 

Earlier approaches on deriving dynamical reductions on stochastic center- 
like manifolds by series expansions are considered by Knoblock and Wiesen- 
feld [IB], Schoner and Haken [27], and Xu and Roberts [33]. 

For stochastic dynamical systems, there are various concepts for invari- 
ance in the definition of invariant manifolds [3J . In the framework of cocy- 
cles [3], the suitable concept for invariance of a random set is that each orbit 
starting inside it stays inside it sample-wise, modulo the change of sample 
due to noise. Another concept is almost sure invariance of a deterministic 
set under stochastic dynamics, i.e., each orbit starting inside it stays inside 
it almost surely. 

In this paper, we consider invariant manifold reductions or restrictions for 
Stratonovich and Ito stochastic differential equations in Euclidean spaces. 

We first study the system of Stratonovich stochastic differential equations 
in R n : 

dX =[AX + F e (X)}dt + B(X)odW(t), X(0)=x , (1) 

where X = X(t,u) is the unknown variable; A is a n x n matrix with k 
eigenvalues of zero real parts and n — k eigenvalues of negative real parts; 
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F e : K — > R and I? : IR n — > M nx?1 are nonlinear vector and matrix functions 
(with e > a small parameter), respectively; and W(t) is a standard vector 
Brownian motion (or Wiener process) taking values in MJ 1 . Moreover, o 
denotes the stochastic differential in the sense of Stratonovich. 

Then we consider the following stochastic system defined by Ito stochastic 
differential equations in ¥L n : 

dX = F(X)dt + B(X)dW(t), X(0) = x , (2) 

where again F and B are vector and matrix functions in M. n and M nxn , 
respectively. And W(t) are standard vector Brownian motion in M n . 

Note that the Stratonovich stochastic differential B(X) o dW(t) and Ito 
stochastic differential B(X)dW(t) are interpreted through their correspond- 
ing definitions of stochastic integrals [2*3] : 

/ B(X) o dW(t) := mean-square lim V B(X( tj+1 ~ * j )) (W t . +1 - W t X 
Jo Atj-to^r ' 2 

3 

[ B (X) dW (t) := mean-square lim V B(X(t j )){W tj+1 - W tj ). 

JO Atj— »0 , 

Note the difference in the sums: In Stratonovich integral, the integrand is 
evaluated at the midpoint tj+1 2 ~ tj of a subinterval (tj,tj + i), while for Ito 
integral, the integrand is evaluated at the left end point tj. See [22] for the 
discussion about the difference in physical modeling by these two kinds of 
stochastic differential equations. There are also dynamical differences for 
these two type of stochastic equations, even at linear level jH]. 

In this paper, we derive an invariant manifold reduction or restriction 
principle for the above systems of stochastic differential equations. 

For the Stratonovich stochastic system (fljl. we consider random invari- 
ant center manifolds. The dynamical reduction is achieved by investigating 
asymptotic behavior of random center manifolds. 

For the Ito stochastic system (J2J), we study deterministic almost sure 
invariant manifolds, which are not necessarily attracting. We reformulate 
the local invariance condition as invariance equations, i.e., first order par- 
tial differential equations, and then solve these equations by the method of 
characteristics. Although the local invariant manifold is deterministic, the 
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restriction of the original stochastic system on this deterministic local invari- 
ant manifold is still a stochastic system but with reduced dimension. 

This paper is organized as follows: In Section [21 we recall some basic 
concepts for stochastic dynamical systems. We consider random center man- 
ifold reduction in Section EJ Finally, in Section HI we construct deterministic 
invariant manifolds by investigating first order partial differential equations 
via the method of characteristics, and thus obtain dynamical restrictions for 
stochastic dynamic systems. 

2 Stochastic dynamical systems 

In this section we introduce some definitions in stochastic dynamical systems, 
as well as recall some usual notations in probability. 

We consider stochastic systems in the state space R n , with the usual met- 
ric or distance d(x, y) = \jYTj=ii. x j ~ Uj) 2 i norm or length ||x|| = y^™ =1 x'j, 
and the usual scalar product < x,y >= Y^j=i x jVi- All invariant manifolds 
and their sample versions are in this state space. 

Some stochastic processes, such as a Brownian motion, can be described 
by a canonical (deterministic) dynamical system (see [3], Appendix A). A 
standard Brownian motion (or Wiener process) W(t) in R n , with two-sided 
time t G R, is a stochastic process with W(0) = and stationary independent 
increments satisfying W(t) — W(s) ~ A/"(0, \t — s\I). Here I is the n x n 
identity matrix. The Brownian motion can be realized in a canonical sample 
space of continuous paths passing the origin at time 



The convergence concept in this sample space is the uniform convergence on 
bounded and closed time intervals, induced by the following metric 



With this metric, we can define events represented by open balls in Q. For 
example, a ball centered at zero with radius 1 is {u> : p(u, 0) < 1}. We define 
the Borel a— algebra T as the collection of events represented by open balls 
A's, complements of open balls, A c, s, unions and intersections of A's and/or 



Q = C (R,R n ) := {u G C(R,R n ) : u(0) = 0}. 




where II ui — to' 



n 
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A c, s, together with the empty event, the whole event (the sample space fl), 
and all events formed by doing the complements, unions and intersections 
forever in this collection. 

Taking the (incomplete) Borel a— algebra T on Q, together with the cor- 
responding Wiener measure P, we obtain the canonical probability space 
(Q, T, P), also called the Wiener space. This is similar to the game of gam- 
bling with a dice, where the canonical sample space is fldice = {1,2, 3, 4, 5, 6}. 
Moreover, E denotes the mathematical expectation with respect to probabil- 
ity P. 

The canonical driving dynamical system describing the Brownian motion 
is defined as 

9{t) : Q -»• Q, 9(t)u(s) := u(t + s) - u(t), s,teR. 

Then 0(t), also denoted as 9t, is a homeomorphism for each t and (t,u) >—* 
9(t)u is continuous, hence measurable. The Wiener measure P is invariant 
and ergodic under this so-called Wiener shift 9 t . In summary, 9 t satisfies the 
following properties. 

• 9q — id, 

• 9 t 9 s = 9 t+s , for all s, t EM., 

• the map (t, uS) \-* 9 t u is measurable and 9 t F = P for all tel. 

We now introduce an important concept. A filtration is an increasing 
family of information accumulations, called cr-algebras, Tt- For each t, a- 
algebra Tt is a collection of events in sample space Q. One might observe 
the Wiener process W t over time t and use Tt to represent the information 
accumulated up to and including time t. More formally, on (Q, T), a filtration 
is a family of a-algebras T s : > s < t with T s contained in T for each s, and 
T s C T T for s < t. It is also useful to think T as the a-algebra generated by 
infinite union of T s y s, which is contained in Tt- So a filtration is often used 
to represent the change in the set of events that can be measured, through 
gain or loss of information. 

For understanding stochastic differential equations from a dynamical point 
of view, the natural filtration is defined as a two-parameter family of cr- 
algebras generated by increments 

Tl := <T{u){n) -co{r 2 ) :s<T U T2< t), s,tER. 
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This represents the information accumulated from time s up to and including 
time t. This two-parameter filtration allows us to define forward as well as 
backward stochastic integrals, and thus we can solve a stochastic differential 
equation from an initial time forward as well as backward in time 

The solution operator for the stochastic system (JTJ or (J2J) with initial 
condition x(0) = x is denoted as p(t,u,x ). 

The dynamics of the system on the state space IR n , over the driving flow 
9 t is described by a cocycle. A cocycle <p is a mapping: 

^:Kxllxl"^ W l 

which is (i3(M) <E> JF ® £>(IR n ), jF)-measurable such that 

tp(0,u,x) = x£R n , 

for t±, ti G K, u G Q, and x G M n . Then ip, together with the driving 
dynamical system, is called a random dynamical system. Sometimes we also 
use <p(t, to) to denote this system. 

Under very general conditions, the stochastic differential systems (0) and 
(J2J) each generates a random dynamical system in M. n ; see IT7] . 

We recall some concepts in dynamical systems. A manifold M is a set, 
which locally looks like an Euclidean space. Namely, a "patch" of the man- 
ifold M looks like a "patch" in M n . For example, curves, torus and spheres 
in IR 3 are one- and two-dimensional differentiable manifolds, respectively. 
However, a manifold arising from the study of invariant sets for dynamical 
systems in R n , can be very complicated. So we give a formal definition of 
manifolds. For more discussions on differentiable manifolds, see 1^3] • 

Definition 1 (Differentiable manifold and Lipschitz manifold) An 

n-dimensional differentiable manifold M , is a connected metric space with an 
open covering {U a }, i.e, M = \J a U a , such that 

(i) for all a , U a is homeomorphic to the open unit ball in R™, 5 = {i£ 
W L : \x\ < 1}, i.e., for all a there exists a homeomorphism of U a onto B, 
h a '■ U a — > B, and 

(ii) if U a nUp ^ and h a : U a —>■ B , hp : Up — >• B are homeomorphisms, 
then h a (U a fl Up) and hp(U a D Up) are subsets ofW n and the map 

h = h a ohp 1 :hp(u a nUp) ^h a (u a nUp) (3) 



(3 



is differentiable, and for all x G hp(U a R Up), the Jacobian determinant 
det Dh(x) ^ 0. 

J/i/ie map (0) is only Lispchitz continuous, then we call M an n-dimensional 
Lispchitz continuous manifold. 

Recall that a homeomorphism of A to B is a continuous one-to-one map 
of A onto B, h : A — > B, such that h' 1 : B ^ A is continues. 

Just as invariant sets are important building blocks for deterministic dy- 
namical systems, invariant sets are basic geometric objects to help under- 
stand stochastic dynamics. Here we present two different concepts about 
invariant sets for stochastic systems: random invariant sets and almost sure 
invariant sets. 

Definition 2 (Random invariant set) A random set M(oj) is called an 
invariant set for a random dynamical system ip if 

ip(t,u),M{w)) C M(9 t cu), teR and oo en. 

Definition 3 (Random invariant manifold) If a random invariant set 
M can be represented by a graph of a Lipschitz mapping 

7*(a>, •) : H + — > H~, with direct sum decomposition H + © H~ = M n 

such that 

M(u) = {x + + 7*(w, x + ),x + G H + }, 
then M is called a Lipschitz continuous invariant manifold. 

We will also consider deterministic invariant sets or manifolds, while the 
invariance is in the sense of almost-sure (a.s.) |H H3J EH I2H1 Ell- 

Definition 4 (Almost sure invariant set) A (deterministic) set M in 
M. n is called locally almost surely invariant for if for all (t ,xo) G M + x 
M, there exists a continuous local weak solution X^ to ' x °^ with lifetime r = 
r(to, xo), such that 

x (to,xa) e Mj V£ > a ^ 

where t At = min(t, r) . 
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3 Random center manifold reduction 



In this section we study the system of Stratonovich stochastic differential 
equations in E n : 



where A is a n x n matrix with k eigenvalues of zero real parts and n — 
k eigenvalues of negative real parts; F 6 : E — > E and B : E n — »• M nxn 
are Lipschitz continuous with Lipschitz constants and Lb, respectively; 
F e (0) = 5(0) = 0. And e is a small parameter so that F e can be seen 
as a small perturbation, that is, we have L e F — > as e — > 0. Let X c be the 
projection of X G E n on ker A (i.e., projection to the center eigenspace). Note 
that later on we will truncate the nonlinearity so that it has global Lipschitz 
constant. For stochastic systems, truncation may not be always appropriate, 
although sometimes it works fine, such as in considering nonlinear dynamical 
behavior near fixed points [7j. We have the following result. 

Theorem 1 (Random center manifold reduction) 

Given the above assumptions for the Stratonovich stochastic differential 
equation Q). If further assume that the equation generates a dissipative 
random dynamical system (e.g., has an absorbing set). Then for small e, 
the long time behavior of can be described by the following stochastic 
differential equation in E fc 



provided is structurally stable. In this reduced system, F € (X C ) := F e (X c + 
0) and B(X C ) := B(X C + 0) with G R n \ ker A, and W c (t) is the projection 
ofW(t) on ker A 

Remark 1 We say the long time dynamics of the stochastic equation (j3J) 
is described by the stochastic equation (jSJ) if both systems have the same 
limit sets (and possibly also share some other invariant sets). 

Remark 2 The random dynamical system tp(t, u) generated by (J3J) is 
called structurally stable, if for any small perturbation (small in the sense 
of the usual metric in the space of continuous functions) to F e (x) and B(x), 



dX = [AX + F e (X)}dt + B(X) o dW(t), X(0) = x , 



(4) 



dX c (t) = F e (X c (t))dt+B(X c (t)) o dW c {t) 



(5) 



S 



the perturbed random dynamical system $(t,u) is tope-logically equivalent 
to tp(t,u). Namely, there exists a random homeomorphism h(u) so that 
$(t,u) o h{uo) = h(6 t (J) o cp(t,u). 

The proof for this theorem can be obtained by modifying the proof in [2H] 
to the systems in the space M n . 
Let us look at an example. 

Example 1 Consider a system of stochastic differential equations in M 2 : 

dx = —xdt + (xy 2 — x 3 — -x)dt + xodWi(t), 
dy = Oydt + (-2 + x 2 y - y 3 - ^y)dt + y o dW 2 (t), 
where W\ and W 2 are independent scalar Brownian motions. 
Let u = (x, y) T . Then 

du = (Au + F(u))dt + Buo dW(t), 

with 

/ -1 \ , (xy 2 -x 3 -\x \ , n „ 1TtjS fxodW 1 (t)\ 
A = n _ J , F(u) = y Q , 2 3 1 and o dW(t) = „ x . L },{ . 

\ J K ' \-2 + x 2 y-y 3 -\y J w \yo dW 2 (t) J 

In order to apply the Ito's formula, we rewrite this system in the equiva- 
lent Ito's stochastic differential equations (see page 36): 

dx = -xdt + (xy 2 - x 3 )dt + xdWi (t) , 

dy = o y dt+(-2 + x 2 y-y 3 )dt + ydW 2 (t), 

where Wi and W 2 are independent scalar Brownian motions. Let u = (x, y) T , 
then 

du = (Au + F(u))dt + Bu dW(t) 

with 

A = ( "o 1 ) = ( -2 X fx~ 2 y-y 3 ) and Bu = ( I ])' 
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Recall the standard scalar product < U\, u 2 >= Xix 2 + y±y 2 and norm ||u|| = 
\/ x 2 + y 2 in R 2 . Then, we apply the Ito's formula (see page 48) to 
obtain "energy" estimate 

^-^-E||ti|| 2 = E(u,du) +F,]-(du, du) 

£ (XL — 

= E(u,du) + E^(BudW{t),BudW{t)) 

= E(u, A(u) + F(u)) + Trace[Bu ■ [Bu) T ] 

= -x 2 + x 2 y 2 - x 4 - 2y + x 2 y 2 - y 4 + ^(x 2 + y 2 ) 

= -^(x 2 + y 2 ) + (2xV - x 4 - y A ) - y + y 2 - 1 + 1 

= -^(x 2 + 2 / 2 )-(x 2 - 2/ 2 ) 2 + ( ?/ -l) 2 -l 

< — E\\u\\ 2 , 
2 ll ll > 

if y is near the equilibrium point (0,0) (so that < y < 1). Note that here 
E denotes the expectation with respect to probability P. This estimate will 
be used to conclude dissipativity for the (truncated) system. 

The nonlinear terms (xy 2 —x 3 ) and (— 2+x 2 y— y 3 ) can be truncated within 
a disk centered at (0, 0) with radius < e < 1 (making them zero outside the 
disk). The truncated nonlinear terms satisfy desired Lipschitz conditions. 
And, the above "energy" estimate implies the dissipative property for the 
truncated system. 

By Theorem Q near the equilibrium point (0,0) (i.e., taking e is small 
enough), the original two-dimensional system is asymptotically reduced to a 
one-dimensional stochastic dynamical system 

dy=(-2-y 3 )dt + ydW 2 (t). 



4 Invariant manifold restriction 

Now we consider the stochastic system (J2J) defined by Ito stochastic differen- 
tial equations in IR n : 

dX = F(X)dt + B(X)dW(t), X(0)=x , (6) 
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where F and B are vector and matrix functions in R n and M nxn , respectively. 
We also assume that F(-) 6 C 1 (M n ;M n ) and B(-) e C 1 (M n ;M nxn ). 

We are going to derive representations of invariant finite dimensional 
manifolds in terms of A, F and B, by using the tangency conditions for a 
deterministic C 2 manifold M in 1": 

fi(u,x) :=F(u,x) ±Y^[DBi(u,x)]BS(u,x) ET X M, (7) 

3 

B>{u,x) E T X M, j = l,---,n, (8) 

where D represents Jacobian operator and Bj is the j—th column of the 
matrix B. The above tangency conditions are shown to be equivalent to 
almost sure local invariance of manifold M; see Filipovic ([EH) and related 

works [ilia mum ig. 

The almost sure invariance conditions (0)-© for manifold M mean that 
the n + 1 vectors, \x and B 3 , j = 1, • • • , n, are tangent vectors to M. Namely, 
these n + 1 vectors are orthogonal to the normal vectors of manifold M. 

In other words, if the normal vector for M at x is N(x), then the almost 
sure invariance conditions ©-© become the following invariance equations 
for manifold M: For all x G M, 

H(x) ■ N(x) = 0, (9) 
B^x)-N(x) = 0, j = l,...,n, (10) 

where • denotes the scalar product in M. n . 

Invariant manifolds are usually represented as graphs of some functions 
in ~R n . By investigating the above invariance equations (|9*|)-([10[). we may be 
able to find some local invariant manifolds M for the stochastic system (JHJ). 

The goal for this section is to present a method to find some of these local 
invariant manifolds. Although the following result and example are stated 
for a codimension 1 local invariant manifold, the idea extends to other lower 
dimensional local invariant manifolds, as long as the normal vectors N(x) (or 
tangent vectors) may be represented; see tangency conditions P|)-([TU|) above 
and O-flHIJ) below. 

Theorem 2 (Local invariant manifold restriction) 

Let the local invariant manifold M for the stochastic dynamical system (0) 
be represented as a graph defined by the algebraic equation 

M: G(x 1 , •■•,£„)= 0. (11) 
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Then G satisfies a system of first order (deterministic) partial differential 
equations and the local invariant manifold M may be found by solving these 
partial differential equations by the method of characteristics. By restricting 
the original dynamical system (QJ) on this local invariant manifold M, we 
obtain a locally valid, reduced lower dimensional system. 

In fact, the normal vector to this graph or surface is, in terms of partial 
derivatives, VG(x) = (G X1 , • • • , G Xn ). Thus the invariance equations (J9|)-(fTUj) 
are now 

fi(x) ■ VG{x) = 0, (12) 
B j (x)-VG(x) = 0, j = l,---,n, (13) 

This is a system of first order partial differential equations in G. We apply 
the method of characteristics to solve for G, and therefore obtain the invari- 
ant manifold M, represented by a graph in state space M n : G(x\, • • • , x n ) = 0. 

In the rest of this section, we first recall the method of characteristics, and 
then work out an example of finding a local invariant manifold and reduced 
system. 

Method of Characteristics: Consider a first order partial differential 
equation for the unknown scalar function u of n variables x±, ■ ■ ■ , x n 

n 

^ ^ G.j(xi, , X n )ll Xi c(xi, , Xn), (14) 

i=i 

with continuous coefficients eij's and c. 

Note that the solution surface u = u(xi, .., x n , t) in x\... ■ ■ ■ x n u— space has 
normal vectors iV := (u xi , • • • , u Xn , —1). This partial differential equation im- 
plies that the vector V —: (oi, • • • , a n , c) is perpendicular to this normal vec- 
tor and hence must lie in the tangent plane to the graph of z = u(x\, • • • , x n ). 

In other words, ( c) defines a vector field in M n , to which graphs 

of the solutions must be tangent at each point [22] • Surfaces that are tangent 
at each point to a vector field in W 1 are called integral surfaces of the 
vector field. Thus to find a solution of equation (JHj), we should try to find 
integral surfaces. 

How can we construct integral surfaces? We can try using the char- 
acteristics curves that are the integral curves of the vector field. That is, 
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X = (xi(t), • • • , x n (t)) is a characteristic if it satisfies the following system 
of ordinary differential equations: 



a x (xi, ■ ■ ■ 



dx. 




dt 
du 

~dt 



c(x\ , . . , x n ) . 



A smooth union of characteristic curves is an integral surface, there may be 
many integral surfaces. Usually an integral surface is determined by requiring 
it to contain (or pass through) a given initial curve or an n — 1 dimensional 
manifold T: 



This generates an n- dimensional integral manifold M parameterized by (s\, .., s n -i, t). 
The solution u(xx, • • • , x n ) is obtained by solving for (s±, .., s n -i,t) in terms 
of variables (xi, ■ • ■ , x n ). 

Remark: If initial data V is non-characteristic, i.e., it is nowhere tangent 
to the vector field V = (ai, • • • , a n , c), and a±, - ■ ■ ,a n ,c are C 1 (and thus 
locally Lipschitz continuous), then there exists a unique integral surface u = 
) containing T, defined at least locally near T. 

Now applying the above method of characteristics to ()12 |) -()13 |) . we obtain 
a solution G = G(x\, • • • , x n ). However, the local invariant manifold M that 
we look for is represented by the equation 



Therefore, a skill is needed to make sure that the solution G = G(xi, • • • , x n ) 
actually penetrates the plane G = in the X\ ■ ■ ■ x n G— space; see Fig. 1. This 
needs to be achieved by sleeting appropriate initial data T. The invariant 
manifold M we thus obtain is defined at least locally near the initial data T. 

We illustrate the method for finding local invariant manifold and the 
corresponding reduced system by an example. 



Xi = f i (s 1 ,..,s n - 1 ),i = l..n 
u = h(s 1 ,..,s n - 1 ) 



G(xi, ■ ■ ■ 



x n ) = 0. 
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Example 2 

dx 

— = x + xWi + xW 2 , 
at 

-j- = 3x + 2y + (x + y)W 1 + (x + y) W 2 
at 

where W{ and are independent scalar Brownian motions. 

We look for a local invariant manifold M C K 2 . For this illustrative 
example, the associated tangency conditions (J7J) and (jSJ) coincide and thus 
becomes a single invariance condition: 

(x,x + y) T e T X M (15) 

We represent the invariant manifold M by G(x,y) = 0. This surface has 
normal vector (G x , G y ). By noticing that normal vector is orthogonal to the 
tangent surface T X M, we see that the above single invariance condition (|15|) 
becomes a single invariance equation: 

xG x + (x + y)G y = 0. 

We solve this first order partial differential equation with initial curve T 
parameterized as (f(s),g(s),h(s)). The characteristic equations are 

dx 

dt = Xl 
dy 

Tt = x + y > 

dt 

We solve these equations and invoke the initial conditions to find that 

x = f{s)e\ 
V = (f(s)t + g(s))e\ 
G = h{s). 

This is the general solution with respect to the general initial condition 
(zo(s),y (s),G Q (a)) := (f(s),g(s),h(s)). By solving for t, s in terms of x, y, 
we obtain G = G(x, y). 
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We illustrate this by a specific choice of initial curve (f(s),g(s),h(s)). 
Note that, in order to obtain a local invariant manifold G(x,y) = 0, we also 
need to pick initial curve so that G actually takes both positive and negative 
values, and thus the invariant manifold G(x, y) = is defined on some set 
in the xy— plane. In other words, the continues function G(x, y) satisfies 
max{G} * min{G} < locally. 

For example, taking T : (x (s), yo(s), G (s)) = (1, s, s), we then have 

V = (t + s)e\ 
G = s. 

Thus s — - — ln{x) and G(x, y) — - — ln(x). Thus an invariant manifold M 
is G(x, y) = 0, i.e., 

— — ln(x) = 0. 

x 
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